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Abstract 

The Lense-Thirring effect is currently being measured by means of a combination of the orbital residuals of 
the nodes of LAGEOS and LAGEOS II and the perigee of LAGEOS II. The claimed total error should be of 
the order of 20%. The most insidious systematic error is due to the mismodelled even zonal harmonics of the 
geopotential and amounts to 12.9%, according to the full covariance matrix of the EGM96 model up to degree 
I = 20. The role and the importance of the LAGEOS-LAGEOS II Lense-Thirring experiment is investigated. 
Using other suitable combinations with orbital elements of the other existing laser-ranged satellites docs not yield 
significative improvements except, perhaps, for one combination including the nodes of LAGEOS, LAGEOS II 
and Ajisai and the perigee of LAGEOS II. The related systematic error due to the mismodelled even zonal part 
of the geopotential reduces to almost 10.7%, according to the full covariance matrix of EGM96 up to degree 
I = 20. 
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I Introduction 

In regard to General Relativity (GR), the experience shows that designing and effectively 
performing experiments aimed to test it is very difficult and often expensive. This is particularly 
true for its linearized weak-field and slow-motion approximation which is adequate in the Solar 
System. The relativistic effects quite often fall below or lie just at the edge of the experimental 
sensitivity of various techniques which could be used. Moreover, there are lots of biases and 
aliasing effects induced by a host of classical phenomena which pose severe constraints to the 
obtainable accuracy in space-based experiments. In many cases the systematic errors are larger 
than the standard statistical errors resulting from the least-squares procedures with which the 
relativistic effects are commonly extracted from the data. 
So, it is of the utmost importance 

• To design experiments capable of effectively testing some post-Newtonian effects not yet 
directly checked at relatively low expense of time and money 

• To assess as more reliably and clearly as possible the error budget so to single out the 
various aliasing competing forces and to evaluate the impact of their systematic errors 

• To investigate more thoroughly and reanalyze the role played by the already existing or 
proposed experiments in order to correctly evaluate their importance 

1.1 The importance of the present Lense-Thirring LAGEOS exper- 
iment 

Nowadays, a first attempt to measure the Lense-Thirring drag of the orbit of a test body 
[Ciufolini and Wheeler, 1995] in the gravitational field of the Earth, by using the following 
linear combination of orbital residuals of the rates of the node 5Q of LAGEOS and LAGEOS 

II and the perigee 6u of LAGEOS II [Ciufolini, 1996] 

5(l l + ciSn 11 + c 2 5u u ~ /i LT 60.2, (1) 

has been reported in [Ciufolini et al, 1998; Ciufolini, 2002]. In eq.(^) c\ ~ 0.295, c 2 ~ —0.35 are 
coefficients suitably designed in order to cancel out every mismodelled contributions of the first 
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two even (I = 2, 4) zonal (m = 0) harmonics J2 and J4 of the multipolar expansion of the non- 
spherical Earth's gravitational field and are built up with the orbital parameters of LAGEOS 
and LAGEOS II (see section 3 for further details). Finally, /xlt is the parameter in terms of 
which the Lense-Thirring drag is expressed: it is 1 in GR and in Newtonian mechanics. GR 
predicts for eq.(|l|) a linear trend with a slope of 60.2 milliarcseconds per year (mas/y). Eq.(Jl]) 
would allow to obtain a claimed total experimental error of the order of 20% [Ciufolini et ai, 
1998; Ciufolini, 2002]. The important GP-B mission [Everitt et ai, 2001], which is aimed to 
the detection of a complementary gravitomagnetic effect on freely falling gyroscopes and which 
is scheduled to fly at the beginning of 2003, would reach a claimed accuracy level of 1%. It is 
worth noting that the strategy of the combined residuals of LAGEOS and LAGEOS II could 
allow to constraining strongly the hypothesis of a Yukawa-like fifth force [Iorio, 2002a] and to 
measure in the field of the Earth the gravitoelectric relativistic perigee advance [Iorio et ai, 
2002a]. 

What is the real importance of the observable of eq.([TJ)? Is it the only possibility we have 
at present in order to measure the Lense-Thirring drag of the orbit of a test body in the 
terrestrial space environment [Casotto et ai, 1990]? Would it be possible to adopt alternative 
combinations involving orbital elements of other SLR satellites? We will try to answer these 
questions in the following. 

In Table 1 we quote the orbital parameters of the existing SLR geodetic satellites and of the 
proposed LARES. In it a, e, i and n are the semimajor axis, the eccentricity, the inclination 
and the Keplerian mean motion, respectively. 

The paper is organized as follows. In section 2 we analyze the systematic errors affecting 
such measurement. In section 3 we look for alternative combinations including the orbital 
residuals of other existing SLR satellites. Section 4 is devoted to the conclusions. 

2 The systematic errors 

The systematic errors affecting such measurement are of two types 

• Long-periodic aliasing harmonics of gravitational (solid Earth and ocean tides [Iorio, 2001; 
Iorio and Pavlis, 2001; Pavlis and Iorio, 2002]) and non-gravitational (direct solar radia- 



3 



Table 1: Orbital parameters of the existing spherical passive geodetic laser-ranged satellites 
and of LARES. Aj=Ajisai, Stl=Stella, Str=Starlette, WS=WESTPAC-1, E1=ETAL0N-1, 
E2=ETALON-2, L1=LAGE0S, L2=LAGE0S II, LR=LARES. a is in km, i in deg and n in 
s" 1 . 





Aj 


Stl 


Str 


WS 


El 


E2 


LI 


L2 


LR 


a 


7,870 


7,193 


7,331 


7,213 


25,498 


25,498 


12,270 


12,163 


12,270 


e 


0.001 





0.0204 





0.00061 


0.00066 


0.0045 


0.014 


0.04 


i 


50 


98.6 


49.8 


98 


64.9 


65.5 


110 


52.65 


70 


n 


0.0009 


0.001 


0.001 


0.001 


0.00015 


0.00015 


0.00046 


0.00047 


0.00046 



tion pressure, Earth albedo, direct Earth IR radiation pressure, thermal effects [Lucchesi, 
2001; 2002]) origin. If their periods P are shorter than the observational time span T Q b s 
they can be viewed as empirically fitted quantities and can be removed from the signal. If 
P = A;T bs, k — 1, 2, ... they average out. If P > T Q b s , as is the case for the insidious 18.6- 
year lunar tide acting on the node and the perigee, the K 1: 1 = 3, m = 1, p = 1, q = —1 
ocean tide, whose period amounts to 5.07 years, acting on the perigee of LAGEOS II, and 
a harmonic constituent due to the direct solar radiation pressure, whose period amounts 
to 11.6 years, acting on the perigee of LAGEOS II as wellQ, they may alias the relativistic 
trend acting as superimposed linear signals. However, the 18.6-year tide is an even zonal 
perturbation which can be canceled out by designing suitably the observable, as done 
with eq. ([[]). The impact of the other two mismodelled harmonics is reduced by the small 
value of the coefficient C2 with which the perigee of LAGEOS II enters the combination 
ofeq.® 

• Secular aliasing trends of gravitational (static even zonal harmonics of the geopotential 
[Kaula, 1966]) and non-gravitational (some tiny thermal effects [Lucchesi, 2001b]) origin. 
The impact of the mismodelled even zonal harmonics is far the most relevant in the error 
budget. Indeed, the systematic error induced by them cannot be suitably handled 
or removed from the signal. One can only assess as more reliably as possible its impact 
on the measurement of /^lt- 

1 Notice that the period of the 18.6-year tide depends only on the luni-solar variables, while the periods of 
the K\ tide and of the direct solar radiation pressure harmonic depend on the orbital geometry of the satellite. 
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Eq. (|I]) cancels out the first two even zonal harmonics: the effect of the remaining higher 
degree mismodelled harmonics 5Jq, 5Jg, ... amounts to 12.9% according to the covariance 
matrix of EGM96 Earth gravity model [Lemoine et al, 1998] up to degree I = 200. 

3 The search for alternative combinations 

A possible strategy for improving the accuracy of the present-day Lense-Thirring LAGEOS- 
LAGEOS II experiment consists of suitable combinations of the orbital residuals 8fl and 8Cj of 
the rates of the nodes and the perigees of different SLR satellites. Such combinations can be 
written in the form 

N 

J2 tkfi = A G r^gr, (2) 
i=i 

in which the coefficients Cj are, in general, suitably built up with the orbital parameters of the 
satellites entering the combinations, the /, are the residuals 5Q, 5io of the rates of the nodes 
and the perigees of the satellites entering the combination, Agr is the slope, in mas/y, of the 
general relativistic trend of interest and /xqr is the solve-for parameter, to be determined by 
means of usual least-square procedures, which accounts for the general relativistic effect. For 
example, in the case of the Lense-Thirring LAGEOS-LAGEOS II experiment [Ciufolini, 1996] 
Alt = 60.2 mas/y. More precisely, the combinations of eq.(||]) are obtained in the following 
way. The equations for the residuals of the rates of the iV chosen orbital elements are written 
down, so to obtain a non homogeneous algebraic linear system of iV equations in N unknowns. 
They are /xqr and the first N — 1 mismodelled spherical harmonics coefficients 5Ji in terms of 
which the residual rates are expressed [Iorio, 2002c]. The coefficients c* and, consequently, Aqr 
are obtained by solving for fx GR the system of equations. So, the coefficients q are calculated in 
order to cancel out the contributions of the first A — 1 even zonal mismodelled harmonics which 

2 In regard to this point, it should be mentioned that, at present, there is no full consensus in the scientific 
community on the reliability of such estimate which might be rather optimistic [Ries et at, 1998]. Indeed, it 
would not be entirely correct to automatically extend the validity of the covariance matrix of EGM96, which is 
based on a multi-year average that spans the 1970, 1980 and early 1990 decades, to any particular time span like 
that, e.g., of the LAGEOS-LAGEOS II analysis which extends from the middle to the end of the 1990 decade. 
Indeed, there would not be assurance that the errors in the even zonal harmonics of the geopotential during 
the time of the LAGEOS-LAGEOS II experiment remained correlated exactly as in the EGM96 covariance 
matrix, in view of the various secular, seasonal and stochastic variations that we know occur in the terrestrial 
gravitational field and that have been neglected in the EGM96 solution. 
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represent the major source of uncertainty in the Lense-Thirring precessions. The coefficients 
Cj can be either constant f] or depend on the orbital elements of the satellites entering the 
combinations. 

Now we expose how to calculate the systematic error due to the mismodelled even zonal 
harmonics of the geopotential for the combinations involving the residuals of the nodes and the 
perigees of various satellites. 

In general, if we have an observable q which is a function q = q(xj), j = 1,2. ..M of M 
correlated parameters Xj the error in it is given by 



5q 



M 



M 

2 S 

h+k=\ \" X h 



dq\ ( dq 



dxi 



a 



(3) 



in which erj = Cjj and al k = Chk where {Chk} is the square matrix of covariance of the 
parameters Xj. 

In our case the observable q is any residuals' combination 



N 



Q = H c ifi(. x i)i j = 1,2. ..10, 



(4) 



i=i 



where Xj, j = 1,2. ..10 are the even zonal geopotential's coefficients J 2 , J±...Jw Since 

dq ^ dfi . 
^— =l^ Ci ^—i 3 = 1,2. ..10, 

dxj ^ dxj 

by putting eq.([|) in eq.(||) one obtains, in mas/y 



(5) 



5q 



10 / N 

i=i 



10 



a 



\i=l 



' dxi 



2 

°hk 



(6) 



h^k=l \i=l 

The percent error, for a given general relativistic trend and for a given combination, is obtained 
by taking the ratio of eq.(D) to the slope in mas yr _1 of the general relativistic trend for the 
residual combination considered. 

The validity of eq.(U) has been checked by calculating with it and the covariance matrix 
of the EGM96 gravity model up to degree I = 20 the systematic error due to the even zonal 



3 In general, the coefficient of the first orbital element entering a given combination is equal to 1, as for the 
combinations in [Ciufolini, 1996; Iorio, 2002b; lorio et ai, 2002a; 2002b]. 
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harmonics of the geopotential of the Lense-Thirring LAGEOS-LAGEOS II experiment; indeed 
the result 

Sfiur = 13% /i L T (7) 

claimed in [Ciufolini et ai, 1998] has been obtained againQ. For the systematic errors due 
to the even zonal harmonics of the geopotential of alternative proposed gravitomagnetic and 
gravitoelectric experiments, according to EGM96, see [Iorio, 2002b; Iorio et ai, 2002a; 2002b]. 
It is worth noticing that, since the orbits of the LAGEOS satellites are almost insensitive to 
the geopotential's terms of degree higher than I = 20, the estimates based on the covariance 
matrix of the EGM96 gravity model up to degree I = 20 should be considered rather reliable 
although the higher degree terms of EGM96 might be determined with a low accuracy. 

It should also be pointed out that the evaluations of the systematic error due to geopotential 
based on the approach of the combined orbital residuals should be free from some uncertainties 
due to possible secular, seasonal and stochastic effects. Indeed, even putting aside the fact 
that most of the secular and seasonal variations of the geopotential are concentrated just in 
the first even zonal harmonics, if we cancel out as many even zonal harmonics as possible, 
the uncertainties in the evaluation of the error based on the remaining correlated even zonal 
harmonics of higher degree should be greatly reduced irrespectively of the chosen time span. 

A very important point to stress is that the forthcoming new data on the Earth's gravita- 
tional field by CHAMP, which has been launched in July 2000, and GRACE, which has been 
launched in March 2002, will have a great impact on the reduction of the systematic error due 
to the mismodelled part of geopotential and on the reliability of the estimates based on it. 

3.1 The combinations without the perigee of LAGEOS II 

In regard to alternative combinations, we first explore the possibility of reducing by 
adopting different residual combinations obtained by substituting the perigee of LAGEOS IlQ 

4 It may be interesting to notice that if, with a more conservative approach, we consider only the diagonal 
part of the EGM96 covariance matrix, up to degree I = 20, the systematic error due to the mismodelled even 
zonal harmonics of the geopotential for the Lense-Thirring LAGEOS-LAGEOS II experiment amounts to 46.5% 
[Iorio, 2002c]. 

5 The orbits of the other SLR satellites are far less elliptical than that of LAGEOS II (except for Starlette), 
so that their perigees are not good observables. It turns out that the inclusion of the perigee of Starlette would 
not yield any improvements. 
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with one or more nodes of Ajisai, Starlette, Stella, WESTSPAC-1, ETALON-1 and ETALON-2. 
The advantages of using only the nodes would be 

• No odd zonal geopotential harmonics 

• No odd degree ocean tidal perturbations 

• Non-gravitational perturbations of smaller magnitude and fairly well modelled 

• Atmospheric drag negligible not only for the LAGEOS and the ETALON satellites, but 
also for Starlette and Ajisai orbiting at lower altitudes 

• Smaller rms on the experimental residuals 

The exclusion of the perigee of LAGEOS II leads to two classes of residual combinations. 

• Combinations without the ETALON satellites. They should be discarded because the 
error due to the uncancelled higher degree even zonal harmonics of the geopotential, 
according to EGM96, is itself far larger than 20%, which is the estimated total error 
in the LAGEOS-LAGEOS II Lense-Thirring measurement. It is so because of the fact 
that the satellites orbiting at lower altitudes than LAGEOS and LAGEOS II are more 
sensitive to the higher degree even zonal harmonics of the geopotential. 

• Combinations with the ETALON satellites. Due to their higher altitude with respect 
to the LAGEOS satellites 5/ilt falls, in some cases, below 1%. However, they should 
not be considered because the coefficients of the nodes of the ETALON satellites are, in 
these cases, of the order of 10 3 so that they greatly amplify every perturbations acting 
on such elements. This is particularly true for the solid Earth tesseral K 1 tide. Indeed, 
it turns out that its I — 2, m — 1, p — 1, q = constituent induces on ETALON-1 a 
perturbation with a nominal amplitude of 1,216.57 mas and period of 10,880.8 days while 
for ETALON-2 we have a nominal amplitude of 1,269.35 mas and a period of 11,130.1 
days. Even if the solid Earth tides are known at a 0.5% level [Iorio, 2001], the effect of 
the too large weighing coefficients and the very long periods of such tidal perturbations 
would induce insidious mismodelled linear trends over reasonable observational time spans 
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of few years. However, it must pointed out that a practical useful use of the SLR data of 
the ETALON satellites would be problematic mainly because the Russian satellites are, 
at present, very poorly tracked. Moreover, they wobble and their center of mass is not 
well defined. Last but not least, the accuracy with which they were constructed and the 
prelaunch ground controls are not comparable to those of LAGEOS satellites. 

In conclusion, even if very appealing, the possibility of using residual combinations involving 
only the nodes of the existing geodetic SLR satellites should be rejected. 

3.2 The combinations with the perigee of LAGEOS II 

It turns out that it is possible to obtain admissible results only by involving the perigee of 
LAGEOS II. 

• Combinations with ETALON-1. We obtained seventeen combinations which, in principle, 
could be considered because < 20% for them. However, it turns out that the 
combinations including the perigee of LAGEOS II and the node of ETALON-1 cannot 
represent a genuine improvement with respect to the present day LAGEOS-LAGEOS II 
Lense-Thirring experiment because £^lt > 13% for them. The systematic errors due to 
the time-dependent non-gravitational perturbations would make the total error budget 
unfavorable with respect to the current experiment. 

• Combinations without ETALON-1. We obtained eight combinations involving only the 
LAGEOS satellites and the lower altitude satellites. It turns out that the combination 

5Q 1 + Cl 8tt n + c 2 5ti A i + c 3 5u 11 ~ ^ LT 61.2, (8) 

with Ci = 0.443, c 2 = —0.0275, c 3 = —0.341 is the only one that could lead to a real, 
slight improvement of the error budget^] [Iorio, 2002b]. Indeed 

— The systematic error due to the even zonal harmonics amount to ~ 10.7% since 
eq.(||) cancels out the effect of SJ G as well 

6 It may be interesting to notice that if, with a more conservative approach, we consider the diagonal part 
only of the EGM96 covariance matrix, up to degree / = 20, the systematic error due to the mismodelled even 
zonal harmonics of the geopotential for such alternative combination amounts to 64.2% [Iorio, 2002c]. 
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— The estimates relative to the periodic time-varying perturbations for eq.([T]) would 
be valid also in this case because the coefficients of the node and the perigee of 
LAGEOS II are almost the same 

— The additional perturbations induced by Ajisai are damped since the coefficient of 
its node amounts only to -0.0275 

— It turns out that the tidal perturbations on the node of Ajisai have periods not 
longer than few years and the non-gravitational perturbations are fairly well studied 
[Sengoku et ai, 1995; 1997] 

— Ajisai, contrary to the ETALON satellites, is a well known and accurately tracked 
geodetic satellite 

4 Conclusions 

The combination of eq.(|l]) plays a really important role in the context of the efforts aimed to 
the detection of the general relativistic Lense-Thirring drag in the gravitational field of the 
Earth. Indeed, at present, it represents a relatively fast, cheap, reliable and almost unique 
way for measuring this effect in the terrestrial space environment by means of the analysis of 
the orbital data of the existing laser-tracked satellites. This should also justify and encourage 
the efforts aimed to improve the error budget and the knowledge of the aliasing classical forces 
affecting such measurement. However, according to the covariance matrix of the EGM96 model 
up to degree / = 20, the alternative combination of eq.(§) might enforce and slightly improve 
the experimental accuracy of the measurement. 

At the present-day level of knowledge of the perturbing forces of the terrestrial space en- 
vironment, only the not yet approved LAGEOS-LARES mission, whose originally proposed 
configuration has been recently modified [Iorio et ai, 2002b], could yield a notable, genuine 
improvement in measuring the Lense-Thirring effect with SLR. 

Moreover, within few years the new data on the geopotential from the CHAMP and GRACE 
missions should be available. They should improve our knowledge, among other things, of the 
spherical harmonics coefficients of the geopotential, and, consequently, the systematic gravita- 
tional errors in all the proposed Lense-Thirring measurements should reduce as well. 
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